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Abstract 

Below the deconfinement phase transition large N c QCD is expected to be a very viscous hadronic 
fluid because both the shear and bulk viscosity to entropy density ratio, n/s,(/s ~ iVf. In this 
letter I show that n/s ~ N% in the confined phase of holographic models of QCD at large N c 
defined in the supergravity approximation. Our results show that the gauge/string duality can be 
used to describe not only nearly perfect fluids but also extremely viscous systems such as a cold 
gas of glueballs. 
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Lattice calculations [1] have shown that the thermodynamic properties of SU(N C ) Yang- 
Mills (YM) theories do not change significantly above T c if N c > 3. The deconfinement 
transition remains first-order and all the thermodynamic quantities, when properly rescaled 
by N%, seem to display a universal behavior when N c > 3. This supports the idea that the 
large N c approximation is a good starting point in the study of the equilibrium properties 
of YM theories and perhaps also the RHIC plasma. It is also natural to assume that the 
large N c approximation can be useful to compute the transport properties of non-Abelian 
gauge theories. 

Indeed, strongly-coupled gauge theories at large N c constructed within the Anti-de Sit- 
ter / Confor mal Field Theory (AdS/CFT) correspondence [2|-4| seem to already capture some 
important aspects of the RHIC data such as the large heavy quark energy loss [5j and the 
nearly perfect fluidity evtdenced by the strong elliptic flow g which — that the de- 
confined medium is some type of strongly-coupled liquid (the sQGP) [7|] where, on average, 
An rj/s ~ 1 Note that after hadronization the matter is expected to be considerably 
more viscous at low temperatures 0, lo| although near T c effects from Hagedorn resonances 
may be able to lower rj/s to a value similar to the one in the sQGP 111 ]. 

It has been shown that rj/s — 1/(Atc) in the deconfined phase of any gauge theory dual to 
supergravity [3] . Even though the deconfined degrees of freedom in these strongly-coupled 
gauge theories do not admit a quasiparticle interpretation, the general idea from kinetic 
theory that strong coupling leads to small mean free path is somehow valid in this case. 
If higher-order derivative corrections are included in the gravity action rj/s changes but it 
remains ~ 1 [3] (from now on ~ 1 means that the quantity does not depend on iV c to 
eading order). At high temperatures, perturbative QCD calculations predict that rj/s ~ 1 



However, below the deconfinement phase transition colorless states are the relevant de- 
grees of freedom and one expects that rj/s ~ N% [k]]. In fact, in QCD at very low 
temperatures momentum isotropization is done mostly by pions and in the chiral limit 
rj/s ~ fn/T 4 ~ N% [15], where f n ~ y/N^ is the pion decay constant. Furthermore, the bulk 
viscosity of a gas of massive pions shows a similar behavior (/s ~ N% 16]. 



In the confined phase of YM at large N c , glueballs are expected to be well defined (non 



relativistic) quasiparticles with negligib 



e width that interact very weakly 



17J. Under such 



conditions, standard kinetic theory [I8j predicts that rj ~ (p) ti£mfp, where the typical 



2 



momentum is (p) ~ \/MT (here M > T is the glueball mass), n is the density, and 
d-MFP ~ A 7 ,? i s ^e mean free path. Indeed, at large N c a gas of glueballs is so dilute and 
weakly interacting that its properties at (and near) equilibrium should be perfectly described 
by kinetic theory. 

While it is by now clear that the deconfined phase of strongly-coupled gauge theories 
with gravity duals exhibits nearly perfect fluid behavior, it is not at all evident that strongly 
interacting models based on holography can also describe the extremely viscous and dilute 
properties of the confined phase of gauge theories at large N c . In this letter I argue that the 
effective colorless degrees of freedom in the confined phase of holographic models of QCD 
at large N c become extremely weakly interacting at low T by showing that r)/s ~ AT?- This 
indicates that holographic methods can also be used to understand the near-equilibrium 
properties of systems that are usually described within kinetic theory. 

In YM at large N c one can expand rj and s as 77 = r/ 2 Af? + r] + O (I/AT?), s = s 2 + so + 
O (1/Af?), where the coefficients rji, Sj depend on T and also the t'Hooft coupling A = g\ M N c . 
In the confined phase s 2 (T < T c ) = while r/ 2 (T < T c ) is nonzero and, thus, rj/s — 772/so 
plus subleading corrections that vanish at large N c . When T/Mq <C 1 (with Mq being the 
mass of the lightest glueball) one obtains that 



M (M T 



3/2 



So(T) ~ ¥ (rirj e ~ Mo/T - (1) 

This approximation gives the smallest so and hence the largest value of rj/s in the confined 
phase. 

At strong coupling, the deconfinement phase transition is understood as a Hawking-Page 
transition [19| from a geometry without a horizon (thermal AdS — > confined phase) to 



another one where a horizon exists (black brane — > deconfined phase) [20(. The transition 
is then first-order with a specific heat of order A,?. Here we assume that the background 
geometry corresponding to the confined phase is gapped and s at low T is given by Eq. 
([1]). Note that while sq in the confined phase is obtained via a one-loop correction in the 
gravity dual, the leading contribution to 77 ~ r? 2 N% at large A is fully determined by classical 
supergravity. 

Some well-known holographic models for the confined phase include bottom-up models 



such as the hard 



2l| and soft wall 



the Sakai-Sugimoto model 



22| models and also top-down constructions such as 



23j. The deconfinement phase transition in these models is a 



Hawking-Page transition [24J, |25| and the confined phase displays a mass gap. However, in 
the hard wall and Sakai-Sugimoto models the masses of very excited hadrons do not follow 
the observed linear Regge behavior M 2 ~ n, where n 3> 1 is the radial excitation number. 



While in the soft wall model [22j linear confinement is built in, the background metric in 
the string frame is identical to AdS§ and, thus, the Wilson loop does not obey an area law. 

Here we will consider holographic models of the type proposed by Gursoy and collabora- 
tors [26] who showed that several properties of YM theories at large N c can be obtained via 
the 5d gravity action (in the Einstein frame) 
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A = [ d 5 xV^G 

167rG 5 J 



n--(d$) 2 - v($) 
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(2) 



which describes the interactions between gravity and a dilaton-like field <£>. The t'Hooft 
coupling in the gauge theory is assumed here to be A = g\ M N c = e*. The basic assumption 
behind such models is that YM at large N c should be described by a noncritical 5d string 
theory where the bulk fields G^ v and $ are dual to and TrF 2 in the gauge theory. The 
axion in the bulk is dual to TrF A F but its contribution to the action is subleading at large 
N r . The effects of the other fields in the bulk are assumed to be somehow included in the 



dilaton potential V(<&). The potential used in [26j is such that the geometry approaches 
AdSs, logarithmically in order to mimic effects from asymptotic freedom. 

It was shown in [26j that a first-order phase transition that closely resembles the decon- 
fining transition in YM at large N c can be obtained if in the IR (where $ 3> 1) the potential 
becomes 

V{<&) ~ -e^ /3 V$. (3) 
Below the critical temperature T c the background fields are the vacuum solutions 

ds 2 = e 2Ao(r) (-dt 2 + dx 2 + dr 2 ) , $ = $ (r) . (4) 

Finite temperature effects are included in Euclidean space t — >■ i r, where r is a periodic 
function r — > t + 1/T. This horizonless geometry corresponds to the thermal AdS solution 
discussed by Hawking and Page in [19( where the boundary is located at r = 0. It is easy to 
see that the leading contribution to the entropy density (~ N 2 ) vanishes in this geometry 
since the system's free energy T = T A does not depend on T. A potential with IR behavior 
given by ([3]) leads to a linearly confined theory with glueball masses that are consistent with 
current lattice data 126] . 
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In hard wall-like models confinement is associated with a scale in the bulk, say, r CT , where 



spacetime literally "ends" 



211 ] . This scale also determines the maximum depth reached by 



the dual string between heavy quark probes at the boundary that experience a linear (hence 
confining) potential. A similar type of singularity below T c is also present in the model ([2]) 
when r — > oo (this is in fact a general property of these two- derivative actions 26] ) but this 
singularity can only be probed by highly excited states while classical string worldsheets 
remain far from this region where the curvature diverges. Above T c , the singularity is 
shielded by the black brane horizon. The maximum depth, r CT , probed by a classical string 
is determined dynamically using the classical solutions of the string equations of motion in 
the background 0I|) and the confining string tension is 



a = T s e 2Ao ^X 4 J 3 



(5) 



where A CT = e* ^ and T s = l/(2ir£ 2 ) is the string tension defined in terms of the funda- 
mental string length. 

The deconfined phase involves a black brane 

dr 2 



ds 2 = e 2A ^ 



-f{r)dt z + dx 2 + 



fir) 



$ = $(r) . 



(6) 



with thermodynamic properties 



26j that match lattice calculations 27]. The horizon func 



tion f(r) has a simple root at the horizon's location r^. Similar black brane solutions where 
investigated in 



28 11 but in these studies the dilaton-like field was assumed to be dual to an 
operator of dimension A < 4 in the boundary. These gravity duals are different than those 
considered in 26[ because they have a nontrivial UV fixed point where the geometry becomes 
AdSz, with radius R and the potential assumes the universal form V($) ~ — 12 + m|$ 2 /2, 
where m$ = A (A — 4) < 0. In this case, $ is dual to TrF 2 at a finite scale in the IR 
above which the system is conformal invariant and not asymptotically free. However, the 
UV behavior is not relevant in our discussion and, in fact, it is sufficient for our purposes 
here to assume that the confined phase is of the form (J4]) with a string tension given by (J5]). 
The shear viscosity in the field theory is computed using Kubo's formula 



V 



limlim / d x -^-([T xy (t,x),T xy (0)}) 



lim lim -ImG fl (u, k) 



(7) 
(8) 
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where G R is the retarded correlator of T xy in the field theory. The graviton absorption 
cross section for extremal black branes was computed in [29(]. According to the AdS/CFT 
correspondence, the energy-momentum tensor is the conserved current that couples to the 
metric disturbances G^ v — > G^ u + 5h^ u at the boundary of the bulk spacetime. In general, 



the line element is assumed to have the diagonal form 



-G m {r)dt 2 + G xx (r)dx 2 + G rr (r)dr 2 . 



(9) 



The fluctuation of the metric relevant to rj is given by Sh xy and in the supergravity approx- 
imation the quadratic part of the graviton action becomes 



62tyG 5 J 



(10) 



where x = $hy. In theories with higher-order corrections a similar action for x can be found 
but the gravitational coupling is in general different than G§ 3(3]. The shear viscosity can be 
obtained from the canonical momentum Il(x, r) associated with the solution of the equations 
of motion for x(x, r) 31 ] 

v ,-,„., , , (ii) 



Il(x, r) 



16tt Gt 



G rr d rX (x,r) 



as follows 



77 = lim lim Im I — — — -r- 

u->Or-+o \ux(w, 0,r) 



(12) 



We shall assume in the following that the disturbance does not depend on the spatial 
momentum and also x(o;,r) = xo(^)9uj( r ) where g u {r) obeys the following equation 



G oG rr Gl x d r ) +u 2 Gl 



9u{r) = 



(13) 



which is obtained from the equations of motion for x- It is convenient 32[ to define the new 
coordinate £(r) = f dr j a/ Gqo G rr G xx so then Eq. ( fl3l) becomes 



(14) 



An important property of Eq. f ll2p is that in the limit u — > (with ux and II fixed) the same 
result for r] is found regardless of the radius used in the evaluation of the right-hand side in 
( fl2|) . as shown in 31] . This follows directly from the definition of II and its corresponding 
equations of motion. Thus, assuming that = £(r*) is well defined at a given r* e (0, oo) 
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and G xx (r*) ^ 0, near r* the general solution of (fl4j) is a sum of incoming and outgoing 

waves at r* 

which can be seen as a consequence of time-reversal invariance. However, in general at finite 
temperature one should expect that time-reversal invariance is broken and the retarded 
Green's function should be associated with the incoming wave at r*, which leads to 

V= 77. n • 16 

In the deconfined phase the black brane contribution to the entropy is simply GHxijh) / (4G5) 
and, after taking r* = r^, one immediately recovers the general result that f]/s = 1/(Att) in 



the deconfined phase of all gauge theories dual to supergravity [12 ]. 

We shall now compute 7] in the background (jlj) that is assumed to describe the confined 
phase of YM at large N c . Such geometries do not have a horizon but the general discussion 
presented above still applies. In fact, rj can also be computed at any radius r* where £* is 
well defined and e j4o<r *- 1 7^ 0. Naively, one may think that in this case a knowledge of the 
full solution of the equations of motion for A (r) is necessary to evaluate rj. However, this 
is not true. In the confining geometries obtained as solutions of (T5]) e A °^ only vanishes at 



the singularity r — > 00 where A (r) ~ — r 2 26) . Moreover, evaluating rj at r a and using ([5]) 
one obtains 

1 la\ m 1 



' 1{T * Tc) = mal {%) at < 17 > 

This is the main result of this letter. It shows that 77 ~ iV 2 in holographic models of YM 
at large N c , which indicates that these models may indeed provide the correct degrees of 
freedom below T c . One should also expect that ( ~ in such geometries. 

Other transport coefficients associated with shear dissipation can be found in the hydro- 
dynamic expansion and to 2nd order in the gradients one finds r n , which basically gives 
the amount of time necessary for a flow gradient to be converted into heat. Causality re- 



quires that (1] / s) / (t^ T) is finite and smaller than unity 33j. Therefore, in order to preserve 
causality in the deconfined phase T n T must be at least of order one while in the confined 
phase it should scale with ~ or even a higher power. Indeed, for a very dilute gas with 



pressure p one can show that T n = rj/p 34] 



It is generally believed that rj/s has a minimum near T c [10(] and, because of the causality 
constraint, one should expect that r n T should also display a minimum near T c . The fact 



that in the confined phase at large N c both r]/s ~ N% and r n T ~ N% or higher indicates 
that the system is so dilute that, in the context of heavy ion collisions, kinetic freeze-out 
would basically occur exactly when the system hadronizes. In this limit, sQGP observables 
such as the elliptic flow coefficient become identically zero below T c . 

One can use Eqs. ([TJ and (1T7|) to derive the upper bound r]/s < rj(T < T c )/s (T) for 
the confined phase of gravity duals. Thus, since this ratio becomes very large at large 
N c we see that the gauge/string duality can be used to describe not only nearly perfect 
fluids but also extremely viscous systems such as a gas of glueballs which, as argued above, 
should be perfectly described using a quasi-particle description derived within kinetic theory 
The results discussed here should also be valid in the case of a more complex holographic 
theory where the confined phase involves not only glueballs but also mesons and hadrons. 
In this case, holographic methods could provide a reasonable framework to investigate the 
properties of the hot hadronic matter produced in relativistic heavy ion collisions. 
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Science Grant No. DE-FG02-93ER40764 and the Helmholtz International Center for FAIR 
(LOEWE program). 
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